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Abstract 

From the literature one infers that the bulk of the order a s corrections to the Drell- 
Yan cross sections da /dm and m 3 d 2 a y 'dmdxp is constituted by the soft plus virtual 
gluon part of the coefficient function. In the case of da /dm it can be shown that at 
fixed target energies the effect of the exact order a 2 corrected coefficient function is 
very well approximated by its soft plus virtual gluon part. Since the complete order 
a 2 contribution to the coefficient function is missing we have to assume that the 
same approximation also holds for m 3 d 2 a/dmdxp. It appears that the discrepancy 
between the exact order ot s corrected cross section and the massive lepton pair data 
taken at fixed target experiments can be partially explained by including the order 
a 2 soft plus virtual gluon part of the coefficient function. 



1 Introduction 



Massive lepton pair production in hadronic interactions is besides deep inelastic 
lepton-hadron scattering one of the most important probes of the structure of hadrons. 
It is well established that one of the dominant production mechanisms is the Drell-Yan 
(DY) process U . Here the lepton pair is the decay product of one of the electroweak 
vector bosons of the standard model (7*, W and Z) which in the Born approximation 
are produced by the annihilation of quarks and anti-quarks coming from the colliding 
hadrons. This process is of experimental interest because it provides us with an al- 
ternative way to measure the parton densities of the proton and neutron which have 
been very accurately determined by the deep inelastic lepton hadron experiments. 
Moreover it enables us to measure the parton densities of unstable hadrons like pions 
and kaons which is impossible in deep inelastic lepton-hadron scattering. Besides the 
measurement of the parton densities there are other important tests of perturbative 
quantum chromo dynamics (QCD) which can be carried out by studying the DY pro- 
cess. Here we want to mention the scale evolution of the parton densities, although 
not observed in this process because of the low statistics, and the measurement of 
the running coupling constant a s (/i 2 ) which includes the QCD scale A. Finally this 
process constitutes an important background for other production mechanisms of 
lepton pairs. Examples are J/^> and T decays or thermal emission of lepton pairs in 
heavy-ion collisions Q. 

The DY process is also of theoretical interest. Since it is one of the few reactions which 
can be calculated up to second order in perturbation theory it enables us to study the 
origin of large QCD corrections which are mostly due to soft gluon bremsstrahlung 
and virtual gluon contributions. In order to control these corrections in the pertur- 
bation series one has constructed various kinds of resummation techniques mostly 
leading to the exponentiation of the dominant terms 0-0. Another issue is the 
dependence of the physical quantities on the chosen scheme and the choice of scales. 
Since the perturbation series is truncated the theoretical cross section will depend on 
the scheme and the renormalization/factorization scale \x. These dependences can be 
reduced by including higher order terms in the perturbation series. An alternative 
way is to determine \i itself (optimum scale) by using so called improved perturbation 
theory like the principle of minimal sensitivity (PMS) ||, fastest apparent conver- 
gence (FAC) [[J or the Brodsky-Lepage-Mackenzie (BLM) procedure [fT0| . 
The first fixed target experiment on massive lepton pair production was carried out 
by the Columbia-BNL group Later on this process was studied in many other 
experiments which were carried out at increasing energies (for reviews see Hl2|). When 
the statistics of the data was improving one discovered that the cross section could 
not be described by the simple parton model given by S.D. Drell and T.M. Yan in 
0. This was revealed for the first time by the NA3 experiment [|T^] (see also [|H{]) 
where the data show a discrepancy in the normalization between the experimental 
and theoretical cross section. This discrepancy is expressed by a so called i^-factor 
which is defined by the ratio between the experimentally observed cross section and 
its theoretical prediction. The above group and the experiments carried out later on 
|I5| show that this fT-factor ranges between 1.5 and 2.5 and is roughly independent 
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of the type of incoming hadrons. The most generally accepted explanation of this 
i^-factor was provided by perturbative QCD. The calculation of the order a s cor- 
rections |IIJ-[]nj to the DY cross section in |l] show that a considerable part of the 
i^-factor can be attributed to next-to-leading order effects. However the order a s 



corrections do not account for the whole i^-factor. More recent experiments p0|-p3 
still indicate that the ratio between the experimental cross section and the order a s 
corrected theoretical prediction is about 1.4, a number which might be explained by 
including QCD corrections beyond order a s as we will show in this paper. 
As has been mentioned at the beginning the DY process is one of the few processes 
where the order a 2 corrections to the coefficient function are completely known. The 
latter refers to the cross section da /dm only where m denotes the lepton pair in- 
variant mass. This coefficient function has been calculated in the MS |23] as well as 
in the DIS |25| scheme. However in the case of the double differential cross section 
d 2 a jdmdxp (d 2 a/dmdy) one has only calculated the order a 2 part of the coefficient 
function which is due to soft and virtual gluon contributions |2£| because the remain- 
ing part is very complicate to compute. Fortunately as is shown in the literature 
[16|]-[[n| the soft plus virtual gluon corrections dominate the total and differential 
DY cross sections in particular at fixed target energies so that we can restrict to 
them to make reliable predictions. 

An analysis of the higher order corrections to the total DY cross section for W- and 
Z-production at large hadron collider energies has been performed in [f24], B^]. Such 
an analysis is still missing for the DY process at fixed target energies and therefore 
we present it here. In particular we want to show that the discrepancy in the nor- 
malization between the order a s corrected DY cross section and the one measured 
at the fixed target experiments can be partially explained by including the order a 2 
contributions due to soft plus virtual gluon effects. 

This paper is organized as follows. In section 2 we present the expressions for the 
various DY cross sections and give a review of the partonic subprocesses included in 
our analysis. In section 3 the validity of the soft plus virtual gluon approximation 
will be discussed and we make a comparison between the order a 2 corrected cross 
section and the most recent fixed target DY data. In appendices A and B we give 
the coefficient functions for d 2 a/dmdxp (d 2 a/dmdy) corrected up to order a s and 
order a 2 respectively. They are presented for arbitrary renormalization and mass 
factorization scale in the MS- as well as in the DIS-scheme. 
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2 Higher order QCD corrections to d 2 a / dmdxp 
(d 2 cr/dmdy) and da /dm 

Massive lepton pair production in hadron-hadron collisions proceeds through the 
following reaction 



H 1 + H 2 -> V + "X" 

U £ 1+ £ 2 (2-D 

Here Hi and H 2 denote the incoming hadrons and V is one of the vector bosons 
of the standard model (7*,^ or W) which subsequently decays into a lepton pair 
(£1/2) • The symbol "X" denotes any inclusive hadronic final state which is allowed 
by conservation of quantum numbers. Following the QCD improved parton model as 
originally developed in || the double differential DY cross section can be written as 



d 2 a 



dQ 2 dxj 



^oy(Q 2 ,M{>) f dt x [ dt 2 H ij (t 1 ,t 2 ,fi 2 )A ij (t 1 ,t 2 ,x 1 ,x 2 ,Q 2 ,fi 2 ). 

i,j •' Xl Jx2 

(2.2) 



Here Q 2 = m? where m denotes the lepton pair invariant mass. The longitudinal 
momentum fraction xf of the lepton pair and the Bj0rken scaling variable are defined 
by 

x F = x 1 -x 2 = —^, r= — = x 1 x 2 , (2.3) 

where stands for the center of mass energy of the incoming hadrons Hi and H2- 
The quantity oy is the pointlike DY cross section which describes the process 

qi + q 2 ^V^£ 1 + £ 2 , (2.4) 

where q\ and q 2 denote the incoming quark and anti-quark respectively. If we limit 
ourselves to V = 7*,Z then oy gets the form 

Ana 2 [ 2 2Q 2 (Q 2 -M 2 ) 
T ee6q — \Z(Q 2 )\ 2 — ete 9 Cv > tCv >* 



w{Q\M 2 ) 



9Q 



(2.5) 



\Z(Q 2 )\ 2 
with 

Z(Q 2 ) = Q 2 -M 2 Z + iM z Y z . (2.6) 

Here the width of the Z-boson is taken to be energy independent and all fermion 
masses are neglected since they are much smaller than y/Q 2 . The charges of the 
leptons and quarks are given by 

2 1 . . 

e* = -l, e u = -, e d = --. (2.7) 



The vector- and axial-vector coupling constants of the Z-boson to the leptons and 
quarks are equal to 

C v , = -C Af (l -4sin 2 ^ 



Ca,£ — 


1 


2 sin29w 


Ca.u = 


—Ca4 = 


Cy,u = 


C A e (l - 
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'V,l — - U W) 

-C A , t (2.8) 

j sin 2 W ) CV jd = -C A/ (l - | sin2 °w) ■ 

The function in (2T2) stands for the combination of parton densities corresponding 
to the incoming partons i and j = q, q, g). Finally Ay denotes the DY coefficient 
function which is determined by the partonic subprocess 

i+ j ->V + "X" (2.9) 

where " X" now represents any multi partonic final state. Both functions i/y and Ay 
depend in addition to the scaling variables U and X{ also on the renormalization and 
mass factorization scales which are usually put to be equal to //. Besides the cross 
section in ( |2.2| ) one is sometimes also interested in the rapidity distribution of the 
lepton pair. In this case the left hand side in (|2.2|) is replaced by d 2 a/dQ 2 dy where 
y denotes the rapidity defined by (see Q2.3|) ) 

y = l\n— , xi = v^e^, x 2 = ^ t~ y (2.10) 



2 x 2 



1 Xp + \/x% + At 

Mn V„ ■ (-MD 



or 



x F + \Jx 2 F + At 

Furthermore on the right hand side the coefficient function Ay is replaced by its 
analogue corresponding to the cross section d 2 a / dQ 2 dy. 

The coefficient function Ay Q2.2|) can be expanded as a power series in the running 
coupling constant a s {n 2 ) as follows 

A, = |(^)"aS"'. (2.12) 

In lowest order the coefficient function of the differential cross section (|2.2|) is deter- 
mined by the subprocess 

q + q^V. (2.13) 

Here V either stands for the virtual photon 7* or the Z-boson and the coefficient 
function is given by 

Ag? = —^—Sfa - x x ) 5{t 2 - x 2 ). (2.14) 
X\ + x 2 

The order a s corrections to the Born process (|2.13|) denoted by A^ are given by the 
one-loop contributions to (2.13) and the gluon bremsstrahlung process 

q + q^V + g. (2.15) 
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In addition to the process above we have another reaction which instead of a quark 
or anti-quark has a gluon in the initial state 



9 + q(q)^V + q(q). (2.16) 

This reaction contributes to A^. Both contributions and A^ have been cal- 
culated in [H], [L8|, |27[ (DIS-scheme) and in (MS-scheme) and are presented in 
( |A.1|) and ( |A.7|) , ( |A.8|) respectively. A part of the order a 2 corrections to the coeffi- 



cient function corresponding to d 2 a/dQ 2 dxp has also been calculated in ||26|| . These 
corrections originate from the soft plus virtual gluon contributions. They consist 
of the two-loop corrections to process ( |2.13| ) and the one-loop correction to process 
( f2.15|) where the gluon is taken to be soft. Furthermore one has also included the 



bremsstrahlungs process 

q + q^V + g + g (2.17) 
and fermion pair production 

q + q + q + q (2.18) 

where the gluons were taken to be soft and the quark-anti-quark pair in the final 
state of ( |2.18|) has a low invariant mass. 



All above corrections contribute to A^ and can be found in appendix B for arbitrary 
factorization and renormalization scale /i where they are presented in the MS- as well 
as in the DIS-scheme. The hard gluon corrections (|2.17| ) and the other two-to-three 
body processes (see below) are very hard to compute at least for the double differential 
cross sections. Fortunately as has been shown in Jl6|]-|[T9|] the bulk of the order a s 
radiative corrections to the cross sections da/dQ 2 and d 2 a / dQ 2 dxp is constituted by 
the soft plus virtual gluon contributions to A^. Therefore within the experimental 
and theoretical uncertainties one can assume that the order a 2 part of the coefficient 
function A q q which is only due to soft plus virtual gluon contributions is sufficient to 
describe the next-to-next-to-leading order DY cross section at fixed target energies. 
This can be tested for the quantity da/dQ 2 which is defined by 

da r x ~ T , d 2 a da r~\^ tf a 



dQ 2 Jt-i dQ 2 dxF dQ 2 -'|inT dQ 2 dy' 

which can also be written as 
da 



Y.MQ\M 2 ) [I ^H,( Xl , X2 ^ 2 )A, (■^-X) (2.20) 



dQ 2 

where A^ now stands for the coefficient function corresponding to the integrated 
cross section da/dQ 2 . 

Since the exact order a 2 corrections to this coefficient function are completely known 
see (MS-scheme) and |25|] (DIS-scheme) one can now make a comparison between 
the exact DY cross section coming from the complete coefficient function and the 
approximate cross section due to the soft plus virtual gluon part. The full order a 2 
contribution to the DY coefficient function requires besides the calculation of the 
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subprocesses mentioned above the computation of the following two-to-three body 
partonic subprocesses. First we have the bremsstrahlungs correction to fl2.16|) 



g + q{q)^V + q(q)+g (2.21) 



which entails the computation of the one- loop corrections to ( 2.16 ). In addition one 
has to add the subprocesses 

qi + g 2 -> V + qi + q 2 (2.22) 

q(q) + q(q) + q(q) + q(q) (2.23) 



and 



g + g^V + q + q. (2.24) 



Reactions P^T| ),(P^D,( P^3D and (ggg ) contribute to the coefficient functions AjJ, 



Ag|, A^ and A^ respectively. The exact result of the coefficient function calcu- 



yy 



lated up to order a 2 for da/dQ 2 gives an indication about the validity of the soft 
plus virtual gluon approximation of d 2 ajdQ 2 dxp (or d 2 a/dQ 2 dy) for which a com- 
plete order a 2 calculation is still missing. In |25| one has made a detailed analysis 
of this approximation for the total cross section of W- and Z- production which is 
derived from ( |2.20| ) by integrating da/dQ 2 over Q 2 . From this analysis one infers 
that the approximation works quite well in order a s as well as in order a 2 when 
My/S > 0.01 provided the DY coefficient function is computed in the DIS-scheme. 
This implies that in practice one can only apply it to the cross section measured 
at the SppS (y/S = 0.63 TeV). The reason that this happens in the DIS-scheme is 
purely accidental. It originates from the large coefficient of the delta-function 6(1 — x) 
appearing in A q g(x) which is small in the MS-scheme. Apparently the combination of 
the anomalous dimension (Altarelli-Parisi splitting function) and the remaining part 
of the coefficient function is very small in the DIS-scheme. It is expected that the ap- 
proximation will even work better when r = Q 2 / S — > 1, a condition which is satisfied 
by fixed target experiments. In this case the phase space of the multi partonic final 
state in the above reactions will be reduced so that only soft gluons or fermion pairs 
with low invariant mass can be radiated off. Their contributions manifest themselves 
by large logarithms of the type (ln fc (l — x)/(l — x)) + which appear in the coefficient 
function in the DIS- as well as in the MS-scheme. 

Notice that the above analysis holds if the mass factorization scale /i is chosen to be 
fi 2 = Q 2 . Therefore it is not impossible that the above conclusions have to be altered 
when a scale completely different from /z 2 = Q 2 is adopted. 

Finally one has to bear in mind that a complete next-to-next-to-leading order analysis 
cannot be carried out yet because the appropriate parton densities are not available. 
The latter can be attributed to the fact that the three-loop contributions to the 
Altarelli-Parisi splitting functions or the anomalous dimensions have not been calcu- 
lated up to so far. Therefore the analysis of the order a 2 corrected result for da / dQ 2 
has to be considered with caution. This holds even more for the order a 2 corrected 
differential distribution d 2 a/dQ 2 dxF or d 2 a/dQ 2 dy. 
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3 Results 



In this section we start with a discussion of the validity of the soft plus virtual gluon 
(S + V) approximation of the order a 2 correction to d 2 a / dQ 2 dx f ( |2.2| ). This is done 
by making a comparison with the integrated cross section da/dQ 2 (|2.19|) for which 
the coefficient function is completely known up to order a 2 . Then we include this 



approximation in our analysis of the fixed target muon pair data published in |20 
||23|| . In particular we show that this correction partially accounts for the difference 
in the normalization between the data in |20[|-||23|| and the order a s corrected cross 



section calculated in [17, 18, 27, 28 



The calculation of the cross sections da/dQ 2 ( |2.19| ) and d 2 a/dQ 2 dxp ( |2.2j ) will be 
performed in the DIS- as well as in the MS-scheme chosen for the coefficient functions 
as well as for the parton densities. The coefficient functions for da/dQ 2 up to order a 2 
can be found in (MS-scheme) and |f25| (DIS-scheme). The coefficient functions for 
d 2 a / dQ 2 dxp corrected up to order a s are obtained from [17|, [18|, [27|] (DIS-scheme) and 



|28|j (MS-scheme). In order to make this paper self-contained we have also presented 
them in appendix A. The order a 2 contribution as far as the soft plus virtual gluon 
part is concerned has been calculated in |^6| and is presented in both schemes in a 
more amenable form in appendix B. For the next-to- leading order nucleon parton 
densities we have chosen the MRS(D-) set [f^J for which a DIS- (A = 230 MeV) 
and an MS-version (A = 215 MeV) exist. Further we use the two-loop (MS-scheme) 
corrected running coupling constant with the number of light flavors nj = 4 and the 
QCD scale is the same as chosen for the MRS(D-) set. For the pion densities we take 
the leading log parametrization (DOl) in |3(J. Using this set one could only fit the 
old lepton-pair data (for references see |30|]) by allowing an arbitrary normalization 
(or i^-factor) with respect to the leading order theoretical DY cross section. In this 
section it is shown that this factor can be partially explained by including higher 
order QCD corrections. Next-to-leading (NLO) order parton densities for the pion 
exist in and |HJ] but they are only presented in the MS-scheme. Also here one 
has to use an arbitrary A'-factor to fit the data which is smaller than found for the 
leading order process since a part of the normalization is accounted for by the order 
a s corrections. Because of the missing (NLO) parton densities of the pion in the DIS- 



scheme we prefer to use the leading log parametrization in |p0| . Finally we choose 
the factorization scale /i to be equal to the renormalization scale where /i 2 = Q 2 . All 
numerical results in this paper are produced by our Fortran program DIFDY which 
can be obtained on request. 

The plots will be presented at three different fixed target energies given by = 
15.4 ; 21.8 and 38.8 GeV/c. At the first energy i.e. y/S = 15.4GeV/c one has observed 
muon pairs produced in the reactions p + W — > + "X" and tt~ + W — > + 

"X" measured by the E537 group The second experiment is carried out at 

y/S = 21.8 GeV/c by the E615 gl| group where the same lepton pair is measured in 
the reaction 7r~ + W — > /x 4 pT +" X" . Finally we discuss the E772 experiment W2 



at \fS = 38.8 GeV/c where the reaction p + N — > + "A" is studied where N is 

either represented by the isoscalar targets 2 H and C or by W (tungsten) which has a 
large neutron excess. Here we will only make a comparison with the 2 f/-data. In the 
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case of the E537, E615 experiments W is given by Z/A = 0.405 whereas E772 used 
tungsten with Z/A = 0.409. Here Z and A denote the charge and atomic number 
of the nucleus respectively. Finally notice that at the above energies we can safely 
neglect the contributions coming from the Z-boson in (|2.5|) since the virtual photon 
dominates the cross section. 

Let us first start with the discussion of the S + V approximation to the coefficient 
function corresponding to da/dQ 2 . The soft plus virtual gluon part of the coefficient 
function, which only appears in A qg , can be written as 



A 



s+v 

qq 



x, Q 



V) 



6(1 



x 



+ E 



An 



'ln 3 (l 



x) 



3=0 



X 



+ 6(1 



x)b®(Q 2 ^ 2 ) 



(3.1) 



where the logarithms have to be interpreted in the distributional sense (see |Ly|). 
The coefficients op and b® depend on Q 2 and the factorization scale fj 2 . The above 
coefficients can be read off the explicit form of ( |3.1| ) given by eqs. (B.3), (B.8) in p4 
and (A. 3), (A. 8) in f25fl . In order to test the S + V approximation to the DY cross 



section we study the following ratios 



R {1) (V^) 



da® da s+v 'W 
+ 



dm 



dm 



da® daW 



(3.2) 



dm 



dm 



and 



daW da® da s+v 'W 
+ —, — + 



dm 



dm 



da® daW da® 
+ — — + 



(3.3) 



dm dm 



dm 



In the above expressions da® /dm (m = y/Q 2 ') denotes the 0(a\) contribution to the 
DY cross section containing the exact 0(a l s ) part of the coefficient function where all 
partonic subprocesses are included. The quantities da s+v ^ / dm stand for the 0(a t s ) 
contribution to the cross sections where only the soft plus virtual gluon part of the 
coefficient function according to (|3.1| ) is taken into account. 

In fig. |l| we have plotted R^(y/r) and R^(y/r) in the DIS-scheme for the \/r-ranges 
explored by the three experiments mentioned above. From the figure we infer that the 
S + V approximation overestimates the exact cross section by less than 10% at small 
^/r-values. At large y^-values this becomes better which is to be expected since in 
the limit r — ► 1 the approximation becomes equal to the exact correction. In this limit 
hard gluon radiation and all other partonic subprocesses like quark-gluon scattering 
are suppressed because of the reduction in phase space. By comparing R®(y/r) 
with R^[y/r) we observe a slight improvement when higher order corrections are 
included in the denominator as well as in the numerator. In fig. |2] we did the same as 
in fig. [l] but now for the MS-scheme. Here we observe that the S + V approximation 
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underestimates the exact DY cross section by more than 10% in particular when 
the CM. energy v^S* is small like in the case of E537 (y/S = 15.4GeV/c) or E615 
(y/S = 21.8GeV/c). Furthermore R^(y/r) Q becomes worse than R [l \^) Q 
in particular in the low y^-region. Hence we can conclude that for da /dm the S + V 
approximation works better in the DIS-scheme than in the MS-scheme. 



In the case of the double differential cross section d a/dmdxF 



m 



y/Q 2 ') the exact 



order a 2 contribution to the coefficient function is not known so that one can only 
make a comparison on the order a s level. The S + V part of the coefficient function, 
of which the explicit form is given up to order a 2 in appendices A and B, becomes 



A s q + v {t 1 ,t 2 ,x 1 ,x 2 ,Q 2 , f x 2 ) 



Xi + x 2 



5(t 1 - xi) 5(t 2 - x 2 ) 




l ki 



(Q 2 ,/i 2 ) 



'ln fc (^/xi - 1)\ (\n\t 2 /x 2 - 1) 



X\ 



x 2 



+S(t 1 -x l ) £ &W(Q 2 ,/i 

k.l 

k+l<2i-l 



2 , / \n k (t 2 /x 2 - 1 



t 2 - x 2 




\ t\ — Xi J 



k.l 

k+l<2i-l 



t l-Xi 



■I'l 



In' 



, 1 -x 2 



x 2 



+5(ti - xi)5(t 2 - x 2 ) 



k.l 

k+K2i 



k 1 Zl± i n ' l ~ X2 



x x 



x 2 



(3.4) 



where the definitions for the distributions indicated by a plus sign can be found in 
appendix A. 

To study the S + V approximation we define an analogous quantity as given for 
da I dm in ( |3.2|) . In the subsequent figures we plot the ratio 



dR^\y/T,x F ) 



d 2 a^ d 2 a s+v ^ 
+ 



dmdxi 



dmdxi 



d 2 a^ d 2 a^ 



(3.5) 



+ 



dmdx p dmdxp 



where the meaning of d 2 a^ / dmdxp and d 2 a s+v,( ^ > / dmdxp is the same as for da^ / dm 
and da s+v '^ / dm defined below ( |3.3| ). Notice that here we cannot present dR^ 2 \\/r, xf) 
because the exact cross section d 2 a^/dmdxF is still unknown. 

Starting with the DIS-scheme we have plotted dR^^y/r^xp) at = 15.4GeV/c 
(E537) for three representative yfr- values as a function of xf in fig. |3|. From this 
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figure one infers that at small yfr around xp = the approximate cross section over- 
estimates the exact one by about 20%. This value is much larger than in the case of 
the integrated cross section da /dm where it was at maximum 10%. The approxima- 
tion becomes better when either \xp\ or yfr gets larger. 

The overestimation is even bigger when the energy increases. This can be observed 
in fig. | (VS = 21.8GeV/c, E615) or fig. | (\fg = 38.8GeV/c, E772). Here one 
overestimates the exact cross section at small ^-values even by 25%. If we repeat 
our calculations in the MS-scheme we observe a considerable improvement of the 
S + V approximation to the double differential cross section (see figs. ||-||). Although 
like in the case of da /dm the approximation underestimates the cross section at high 
yfr- values the difference with the exact one is less than 5%. 

Summarizing our findings we conclude that in the case of the DIS-scheme the S+V ap- 
proximation works better for da /dm than for d 2 a / dmdxp whereas for the MS-scheme 
just the opposite is happening, except for r — > 1 where R^(y/r) and d R^\y/r, Xp) 
become close to 1 independent of the chosen scheme. Further from figs. |3|-|]it appears 
that when d 2 a s+v '^ /dmdxp is integrated over xp according to ( |2.19| ) we get a result 
which differs from the one obtained from da s+v '^ / dm in (|2.20| ) in particular at small 
yfr. On the first sight this is surprising because one expects the same cross section 
da /dm independent of the order of integration. However both procedures only lead 
to the same answer for da /dm when the full coefficient functions are inserted in the 
equations for d 2 a / dmdxp ( |2.2| ) and da /dm ( |2.20| ). If we limit ourselves to the S + V 
part of the coefficient functions as given in ( |3.1|) and (|3.4j) then the two procedures 
to compute da /dm only provides us with the same answer when r — > 1. This we 
have also checked for the order a 2 S + V contribution. Therefore the expression in 
(|3.1|) is not the integrated form of equation ( |3.4| ) except if r — > 1. This explains why 
at large r i? (1) (v^) O and dR^(y/r, xp) ( p.5| ) are roughly the same and equal to 
1 irrespective of the chosen scheme. The above properties of the S + V approxima- 
tion also reveal that if yfr becomes much smaller than 1 one has to be cautious in 
predicting the still unknown d R ( - 2 \y/r , x p) from the values obtained for the known 
Rw(y/r) ( j3.3|) and d R^(y/r, xp) (|3.5|). In the subsequent part of this work we will 
use as a guiding principle that as long as \d R^^y/r, xp) — 1| < 0.1 we expect that 
the S + V approximation of the second order contribution to d 2 a/dmdxp will be 
very close to the exact result. If \d R^'(\/t, xp) — 1| > 0.2 then one should not trust 
this approximation and one has to rely on the predictions obtained from the first 
order corrected cross section. This implies that for the experiments discussed in this 
paper one can make a reasonable prediction for the second order correction as long 
as y/r > 0.3. 

After having discussed the validity of the above approach at fixed target energies we 
will now make a comparison with the data of the E537 |2(|, E615 [^TJ and E772 |22|, |23 



experiments. For that purpose we compute the Born cross section d 2 ao/dmdxp, the 
order a s corrected exact cross section d 2 a\j ' dmdxp and the order a 2 corrected cross 
section d 2 a2/ 'dmdxp. Notice that in the latter only the contribution due to the co- 
efficient function k- s q £ v ( |3.4j ) (see appendix B) has been included because the other 
contributions are still missing. The computations have been carried out in the DIS- 
scheme. The results for the MS-scheme will be shortly commented upon at the end 
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of this section. 

Starting with the experiment E537 (y/S = 15.4GeV/c) we have plotted the quantity 

~T~ = I dxp T~7~ ( 3 - 6 ) 
dm Jo dmdxp 

in figs. | and |10j for the reactions p + W -> + "X" and ir~ + W -> fi+fi~ + "X" 

respectively. Notice that xp in |20| is defined as xf = 2pz,/[(l — r)y/S] which differs 



from the usual definition in []I7|, 0, [?7|, Since the higher order QCD corrections 
are calculated for d 2 ajdmdxp with xp defined in (|2.3| ) and the cross section is not 
a Lorentz invariant we had to change the rcp-bins in table III of [20| according to 



our definition above. Figs. || and |10| reveal that the data are in agreement with the 
order a s as well as with the order a 2 corrected cross section but lie above the result 
given by the Born approximation. The difference between the latter and the data is 
observed when we consider the quantity 

da r 90 d 2 a 

dm- — - — (3-7) 



dxp J 4.o dmdxp 

which is presented in figs. [II] and |T2] for the above two reactions. Even the order a s 
corrected cross section lies below the data for xf < 0.6 as can be seen in fig. [12|. On 
the other hand the order a 2 corrected cross section is in agreement with experiment 
over the whole xp range. 

The second experiment, E615 |2l] also studies the reaction n~ + W — > fi~fi + + "X" 
but now for y/S = 21.8 GeV/c. In fig. [13] we have compared the quantity da/dy/r = 
\/S da/dm with the data where da/dm is defined in the same way as in (|3lf ). Apart 
from the bump, which is due to the T resonance at about \/r = 0.43, the order 
a 2 corrected cross section reasonably describes the experimental results whereas the 
Born and the order a s prediction fall below the data. The importance of the order 
a 2 contribution is also revealed when we study the double differential cross section 

d 2 a 1 fV^ , d 2 o 

dV^ , 7=^ (3.8) 



d^frdxp y/T2 — \fr\ J yrf d^frdxp 
for various xp regions, see figs. |14|4T9]. The curves predicted by the Born and the order 



ot s corrections all lie below the data. For ^pr > 0.277 even the order a 2 contribution 
is not sufficient to close the gap between theory and experiment. This is due to the 
presence of the T in the region 0.323 < \/t < 0.599 which has not been subtracted 
from the data. The discrepancy between the order a 2 corrected cross section and the 
data becomes even more clear when we plot the i^-factor (fig. ^) defined by 

, d 2 a, 
dxi 



F d^frdxp 



o 
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in fig. ^0] and compare the above expression with the experimental J\~-factor which is 
given by 

dxp 



Aexp(^) = d ~^~ (3.10) 



/ 

JO 



dx 



F 



dyfrdxF 

where d 2 Oij 'dy/rdxF denotes the order a\ corrected cross section. Fig. EU] shows that 



neither K\ nor K2 fit the data. The second order corrected J\~-factor is closer to the 
data in the small v^r-region. It is a pity that due to the presence of the T in the 
data it is difficult to compare theory with experiment in particular in those regions 
of y/r where the S + V approximation is supposed to work. 

Finally we also made a comparison with the data obtained by the E772 experiment 
for the reaction p + 2 H — > + "X" carried out at v^S* = 38.8 GeV/c. The main 

goal of this experiment was to find a charge asymmetry in the sea-quark densities of 
the nucleon i.e. u(x) 7^ d(x). Here we are also interested whether the data obtained 
for m z d 2 a / dmdx f are in agreement with the order a 2 s corrected DY cross section. In 
fig. |2l] we have plotted the data for m = 8.15 GeV/c 2 and compared them with the 
predictions given by the Born, the order a s corrected and the order a 2 corrected cross 
section. The figure shows that the order a 2 corrections are needed to bring theory 
into agreement with the data. Notice that at this m-value one obtains \fr = 0.21 
which is quite small for the S + V approximation so that the result has to be in- 
terpreted with care. In the next figure (fig. ^) we study the effect of the higher 
order QCD corrections on the suppression of the cross section near xp = 0.0 which 
is caused by the difference between the up-sea and down-sea quark densities. Notice 
that the pp reaction is symmetric whereas the pn reaction is asymmetric around 
xp — 0.0 irrespective whether there is charge asymmetry or not. Therefore the pn 
reaction leads to an asymmetry even for isoscalar targets like 2 H . In fig. ^ we 
have presented the order a 2 corrected cross section for three different parton density 
sets for the nucleon. They are given by MRS (SO) and MRS (DO) where the former 
has a symmetric sea (u(x) = d(x)) whereas the latter contains an asymmetric sea 
(u(x) 7^ d(x)) parametrization. For comparison we have also shown MRS(D-) which 
only differs from MRS (DO) that the gluon and sea densities have a much steeper 
small x-behavior (lipatov-pomeron) than the ones given by MRS (DO) and MRS (SO) 
(non perturbative pomeron). Fig. [22| reveals that there is hardly any suppression of 
the cross section for xf < while going from the symmetric sea (MRS (SO)) to the 
asymmetric sea (MRS (DO)) parametrization so that both parton density sets are in 
agreement with the data. 

If other parton densities are used like those discussed in [Q the suppression for 
Xf < can be much larger. For the MRS-set it appears that a change in the small 
x-behavior of the parton densities leads to a larger suppression of the cross section 
(compare MRS (DO) with MRS(D-)) than the introduction of a charge asymmetry in 
the sea-quarks (MRS(SO) versus MRS(DO)). 

In addition to the calculations performed in the DIS-scheme we have also presented 
in figs. 0-^l| the order a 2 corrected cross section computed in the MS-scheme. Al- 
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though the latter is an improvement with respect to the order a s corrected result 
it is smaller than the cross section computed in the DIS-scheme except when xp is 
large. This is not surprising because figs. |5|-[8] already indicate that the approximation 
underestimates the exact cross section in the case of the MS-scheme. 
Summarizing the content of this work we can conclude that up to the order a s level 
the soft plus virtual gluon contribution gives a fairly good approximation of the exact 
DY cross section d 2 a j dmdxp. Therefore we expect that this approximation will also 
work for the a 2 correction as long as the cross section is computed at fixed target 
energies and for yfr > 0.3. In this r-region we expect that all other partonic sub- 
processes are suppressed due to the reduction in phase space. This expectation is 
corroborated by a thorough analysis of the second order contribution to da /dm for 
which the exact coefficient function is known. Because of the missing pieces in the 
order a 2 contribution to the coefficient function corresponding to the cross section 
d 2 a / dmdxp and the absence of the next-to-next-to- leading order parton densities we 
have to rely on the order a 2 soft plus virtual gluon approximation to make a compar- 
ison with the data. Using this approach we can show that a part of the discrepancy 
between the data and the order a s corrected cross section can be attributed to the 
higher order soft plus virtual gluon contributions. 
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Appendix A 



In this appendix we will present the order a s contributions to the coefficient functions 
corresponding to d 2 a/dQ 2 dxF coming from the partonic subprocesses in ( 2.15| ) and 
( P-16|) . Although these processes have been calculated in the DIS-scheme in |L7|, [L8| 
(see also \QTn) and the MS-scheme [281 we have some different definitions for the 



distributions and we have a small disagreement with the coefficient function for the qg 



subprocess in [g8j . Moreover we want to give a clear definition for the soft plus virtual 
(S + V) gluon part of the coefficient function corresponding to the qq subprocess. 
We have recalculated the double differential cross section d 2 o jdQ 2 dxp for the partonic 
subprocesses ( |2.15| ) and ( |2.16| ). After performing the mass factorization in the MS- 
scheme the coefficients (see the definition in ( |2.12j )) read as follows 



^qq (^1) Xi,X 2 , Q 2 , /i 2 ) — Cj 



Sfa - xi)5(t 2 — x 2 ) 
x± + x 2 



4 In 



(l- Xl )(l-x 2 ) 

X\X 2 



+ 6l In % - 16 + 12C(2) + 2 ln ; 
I fi 2 



(1 - Xi)(l - x 2 ) 



X\X 2 



Xi 



Xi + x 2 



X'2 



/i 2 \h 



x 2 J + 



In 



1 — X\ 

X\ 



+4 



t 2 - X 2 



4 xi + x 2 t 2 
+ In — + 2— 



x 2 



x 2 



X\ + t 2 



t 2 



2 t2 + °° 2 In — 1 +X2 ^ 1 _Xl )( t2 ~ X2 "> 



' 9 



1 



XxX 2 {t 2 + Xx) 

4 



+ [h 



(tl - Xi) + {t 2 - x 2 ) + 

1 \ 4 



X\ + x 2 
h + xi 

t 2 2 (h + x 2 ) ~ t\{t 2 + Xl ] 



- 2 



t 2 + x 2 



+ 



* h,x 1 

1 

t 1 -x 1 J + 
2 



x 2 \ 



+ 



t 2 -x 2 J + (t 2 + Xi)(ti + X 2 ) t 2 (ti + X 2 ) ti(t 2 + Xi) 



2xo 



+ 



2(x 1 +x 2 )(t 2 1 +t 2 2 ) 



(a.i; 



where the color factor C F is given by C F = (N 2 - 1)/2N (QCD : N — 3). In 
this appendix and in the next one the distributions indicated by a plus sign in the 
denominator are defined as 



*k V tk — Xk I , 



1 dt k W{tk/Xk 1} (/(t fc ) - f(x k )) (A.2) 



/■i /-i h\\t 1 /x l - 



'\n j (t 2 /x 2 



x 2 



f{ti,h 
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dti / dt 2 7 t lf{ti,h)-f(xi,t 2 ). 

xi Jx2 t\ — X\ t 2 — X2 



-f(ti,x 2 ) + f{x 1 ,x 2 )) 



(A.3) 



Expression (|A.l ) for is in agreement with eq. (A. 4) in [28|. Notice that the au- 



thors in |28[ give a different definition for the distributions. This leads to a difference 
between ( |A.2 ) and eq. (A. 12) in which equals 

lnt fe /x fc 



dt fc - 



-In 2 x k + Li 2 (l - x k ) 



where the dilogarithmic function Li 2 (x) is defined by 



U 2 (x) = - -ln(l-t) 
Jo t 



(A.4) 



(A.5) 



The expression between the square brackets in (|A.4j) , multiplied by two, has to be 
added to the coefficient of the S(ti — Xi)5(t 2 — x 2 ) term in eq. (A.4) of [^pj so that 
one obtains the same result as we have in ( |A.1| ) above. 

The soft plus virtual gluon part of is defined by isolating the double singular 
terms in flOJ ) of the types 5(ti-xi)5(t 2 -x 2 ), S(t x -xi) (5=^) , 5(t 2 -x 2 ) (j^) 

and \-r^ — ) (r^ — ) • Hence we obtain 

5(ti - xi)8(t 2 - x 2 ) 



A 



S+V,(l) 
11 



{h,t 2 ,x x ,x 2 ,Q ,/i ) = C 



4 In 



(1 - Xi)(l - x 2 ) 



6 \ In 



+4 



XiX 2 

5(U - xi) 

Xi + x 2 

/ \n(t 2 /x 2 - 
v t 2 — x 2 



Xi + x 2 



16 + 12C(2) + 21n i 



(1 - Xi)(l - x 2 ) 



X\X 2 



1 ®\ 

t 2 — x 2 / + [i \l 2 — x 2 



+ [tl <-> t2,^l ^2] 



hi 



1 — X\ 

Xi 



Xi + x 2 



(ti - Xi)+(t 2 - £2)+ 



where we have taken the residues at t k = x k . 

For the gq subprocess we obtain the coefficient function 



AW{t 1 ,t 2 ,x 1 ,x 2 ,Q 2 ,/i 2 ) = T J 



S(t 2 - x 2 ) 
Xi + x 2 



2 *i + (* 



'1 



ln 



Q 2 
/i 2 



| ^ x 2 + (ti - xi) 2 ^ (xi + x 2 )(l - x 2 )(ti - xi) t 4xi(ti-xi) 



t? 



XiX 2 (ti + x 2 ) 



tl 



(A.6) 
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-4 



XiX 2 (x 1 + x 2 ) 

7372 
L l L 2 



i X l ~ X 2 + 2 1^ + X2 )( t2 ~~ ^2)^2 



tH 



1 L 2 



^1^1(^1 + ^2 



+4 



and 



A^ ) (ti,t 2 ,a;i,a;2,Q 2 ,/i 2 ) = A g 1 q ) (t 2 ,t 1 ,x 2 ,x 1 ,Q 2 , 



(A.7) 



(A.8) 
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where Tf = 1/2. 

There is a discrepancy between our answer in (|A.7|) and the one given in eq. (A.8) of 

x'^ -\~(t-[ X^ ) ^ 

The difference between their result and ours equals 2 1 a — — . This discrep- 



ancy can be attributed to the procedure that in n- dimensional regularization before 
mass factorization the cross section with one gluon in the initial state has to be di- 
vided by n — 2 in order to average over the initial gluon polarizations. Only in this 
case one can combine the coefficient functions with the parton densities of which the 
scale evolution is determined by the two-loop anomalous dimensions (or Altarelli- 
Parisi splitting functions) calculated in the literature (see e.g. fl32"f ). The expression 
in eq. (A.8) of |28| can be only obtained if the polarization average factor is a 1/2 
instead of l/(n — 2). In the latter case one has to modify the two-loop anomalous 



dimensions via a finite renormalization. However the MRS parton densities in |29 



were constructed using the anomalous dimensions in p2j so that one has to divide 
the parton cross section by n — 2 and not by 2. The choice of the polarization aver- 
age factor shows up again when we want to present the coefficient functions in the 
DIS-scheme. The results in the DIS-scheme are obtained by performing a finite mass 
factorization. The coefficient functions in the two schemes are related by 



A-ij(ti, t 2 , Xi,x 2 , Q 2 , fj, 2 



DIS 



E 

kJ 



1 ,1 

dux I du 2 Y ki 

XI Jx2 



U2 



A M (ui,u 2 ,xi,x 2 ,Q 2 ,fi 2 ) 

Up to order a s , T qq (x) and T qg (x) are given by 

On(l — x) 



T qq (x) 



5(l-x) + -^C F 
An 



x 



MS' 



(A.9) 



2(1 + x) ln(l - x) + 6 + Ax 



F qg (x) 



_ 2 i±^lnx- (VM +5(l-x)(-9-4C(2)) 
I — x VI — xj + 



J- f 

An 1 



2 \x 2 + (1 - x) 2 } In + 16x(l 

L ' x 



(A.10) 



(A.ll) 



Expressions ( |A.10| ) and ( |A.11| ) are in agreement with Cp\ an d Cq \ in appendix I of 
|33H . Notice that the authors in f28j used a T qg (x) where 16x(l — x) — 2 is replaced 
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by 12x(l — x) which is obtained when the gluon polarization average factor is taken 
to be 1/2 instead of l/(n — 2). See the discussion above. 
The coefficient functions in the DIS-scheme read 



Atf(ti,h,x 1 ,x 2 ,Q' 2 ,n 2 ) = C F 



{l-x 1 )(l-x 2 ) 



XiX-2 



8{t i - xi)5(t 2 - x 2 ) 

Xi + X 2 

+ 6lln^ + 2 + 20C(2) + 31n 



41nl^lnl^ 



Xi x 2 

(l-x 1 )(l-x 2 ) 



A* 



x x x 2 



5(ti - xi) 
Xi + x 2 



-!41n^i + 3 

Xi 



t 2 - x 2 

1 



t 2 - x 2 



+ 



t 2 + x 2 
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1 



In 



X\ + x 2 
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t 2 



1 



Xi + x 2 
t\ + Xi 



Xi(t 2 +x ± ) 
4 



t 2 - X 2 Xi+ t 2 

+ [*1 ^2, ^1 ^2] 
t 2 + ^2 



£2 i| 



(tl - Zl)+(t 2 - x 2 )_ 



1 



+ 



ti - x 1 
2 



t 2 -X 2 J + (t 2 + + X 2 ) t 2 (ti+X 2 ) tifo+Xi) 



-2- 



- 2- 



£ 2 



+ 2 



(rri +x 2 ){t\ + tl) 



(A.12) 



tl(*i + x 2 ) t?(t 2 + a;i) + 1 2 ) 

The soft plus virtual gluon part is obtained in the same way as discussed in the case 
of the MS-scheme 

5(t l - x l )5(t 2 - x 2 ) 



A^ y ' (1) (ti,t 2 ,x 1 ,x 2 ,Q 2 , Ai 2 ) = C f 



Xi + x 2 



4 In In I^ 2 - 
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Xi x 2 
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X\X 2 
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X\ 
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1 



+C i 



1 



Xi + x 2 



t 2 - x 2 
4 



+ [ti <-> t 2 , Xi <-> x 2 ] 



(ti - ari)+(t 2 - a; 2 )+ 



(A.13) 



The coefficient function for the subprocess with the gluon in the initial state becomes 



A$ {tl ,t 2 ,x 1 ,x 2 ,Q 2 ,v 2 )=T f 6{t \ X2) [{2 

x\ -\- x 2 yy 



x\ + (h-x^ 



In 
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+2 



x\ + (ti - X\ 



In 



(xi + x 2 )(l - x 2 ) 
X 2 (t\ + x 2 ) 



2 _ 



£1 + x 2 



*2 - ^2 



X2 + (tl + X 2 ) 2 Xi + X 2 



tl(t 2 + 2l) 



+ 4- 



-4 



+4 



XlX 2 (Xx + x 2/ 
l l l 2 

XiX 2 (Xi + x 2 ) 



_ | x l x 2 
*f«2 



+ 2 



(Xi + X 2 )(t 2 -^2)^2 + ^l) 



(A.14) 



where A^ is related to A^ via relation ( |A.8|) . 

We have explicitly checked that if the above coefficient functions are inserted in ( p.2| ) 
and the integrals over xp are performed according to ( [2.1 9|) one gets the same answer 
as given by da /dm ( |2.20| ) with the coefficient functions obtained from [21, ^5|. The 
coefficient functions for d 2 a/dQ 2 dy have not been explicitly listed here but are present 
in our computer program D1FDY. In the case of d 2 a / dQ 2 dy we agree with the results 



for the MS-scheme published in [28] except for l/2t\ in eq. (A. 20) which has to 
be replaced by X1 \ tl ~ xx ) , This difference follows again from taking the average over 

the initial gluon polarizations as discussed for A^ above. Our results for the DIS- 



scheme agree with those presented in the appendix of [27|. Notice that the soft plus 



virtual gluon part of A^ for d 2 a/dQ 2 dy can be obtained from (|A.6|) and (|A.13|) by 
multiplication with x\ + x 2 . 
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Appendix B 



The order a 2 contribution to the coefficient function in the S + V approximation 
has been calculated in P5fl . Including the mass factorization parts represented by 
\\iQ 2 / li 2 and rewriting the coefficient function in a more amenable form as presented 
for the first order correction in appendix A it reads in the MS-scheme as follows 
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+n/C F { 



-DoitjDoih) 



L M + -(£>i(ti)A>(* 2 ) + A>(*i)£>i(*2)) 



40 1 

— P> (*l)A)(* 2 )} 



(B.l) 



Here the color factors are given by C A = N,C F = (N 2 - 1)/2N (QCD : N = 3) and 
rif denotes the number of light flavors. In the above expression we have introduced 
the following shorthand notations 



P\x k ) = In 
Di(t k ) = 



1 ~ %k 

x k 



tk %k 



+ 



4 = ln * % 

In the DIS-scheme the above coefficient function becomes 
A s+v,(2) ( . . 2 2 5(t 1 -x 1 )5(t 2 -x 2 ) 

L±qq [tl,t 2 ,Xl,X 2 ,LJ ,1-1) — — 



(B.2) 
(B.3) 
(B.4) 



Xi + x 2 



18 - 8C(2) + 8 (P 2 ^) + P 2 [x 2 )) + 24 (P( Xl ) + P{x 2 )) 



+ 16P(xi)P(x 2 ) 



4 + 



15 + 84C(2) + 48C(3) + 12 (P 2 { Xl ) + P 2 {x 2 )) 



+ 16 (p 2 (xi)P(a; 2 ) + P(xi)P 2 (a; 2 )) + A8P( Xl )P(x 2 ) 



964 



L M + 14C(2) + 72C(3) + ^C 2 (2) 



+ (26 + 64C(2))(P(x 1 ) + P(x 2 )) 
+12 (P 2 ( Xl )P(x 2 ) + P( Xl )P 2 (x 2 )) + 8P 2 (x!)P 2 (x 2 ) 
+ (17 + 80C(2))P(x!)P(x 2 ) + ( 9 - - 8C(2)) (P 2 ^) + P 2 (* 2 )) 
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+ (y + 36C(2) + 24C(3)) (P( Xl ) + P(x 2 ))] 



+CaCf 
22 



22 

-u-—(p(x 1 ) + p(x 2 )) 



-y + ^)) + (— - 8C(2)j (P(xO + P(x 2 )) 



1 93 

— - 22C(2) - 24C(3) 



44 



P(*i)P(* 2 ) 



L> + *f + ^C(2)-36 C (3)-fc 2 (2) 



-y (P 2 (^i)P(^ 2 ) + P(^i)P 2 (^ 2 )) - y (P 2 (^i) + P\x 2 )) 



2 

+ (y - 12C(3)) (P(*i) + P(x 2 )) + ( 2 f~ 8C(2)) P(zi)P(* 2 )} 



+nfCj, 



2 + -(P( Xl )+P(x 2 )) 



-f + 4C(2) 



+ ^ (P 2 ^) + P 2 (X 2 )) + ^P(X!)P(X 2 ) - + P(X 2 )) 



3 

380 
"IT 



P„- 



38 



4D i 

C(2) - -P^OP^a) + - (P 2 { Xl )P{x 2 ) + P^OP 2 ^)) 



+ (P 2 ( Xl ) + P 2 (x 2 )) - 5(P(xi) + P(x 2 )) 
<5(ti - xi) 



Xl + x 2 



C 2 



16pi(t 2 ) + P(x 1 )D (t 2 )) + 2AD (t 2 ) 



+ 



32P(x 1 )D 1 (t 2 ) + 16P 2 (x 1 )D (t 2 ) + 24P>!(t 2 ) + 48P(a; 1 ) J D (t 2 ) 



+ (26 + 64C(2))D (t 2 ^ 



Lfj, + 16P 2 (x 1 )D 1 (t 2 ) + (9 - 16C(2))A(t 2 ) 



+24P(x 1 )D 1 (t 2 ) + 12P 2 (a: 1 )r, (t 2 ) + (17 + 80C(2))P(x 1 )D (t 2 )) 
15 



+ 



— + 36C(2) + 24C(3))r, (t 2 )} 
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-yA)(* 2 ) 



44 44 
-yP>i(* 2 ) - —P(x 1 )D (t 2 ) 



9 



-8C(2) P» (t 2 ) 



44 22 
L M - —P^OACta) - -P 2 ( Xl )D (t 2 ) 
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II AIM ( ^ - 8C(2)) P(xi)D Q (t 2 ) + (y - 12C(3)) A,(* 2 )} 



|a>(* 2 ) 



|(D!(t 2 ) + P(xi)D (t 2 )) -yA)(* 2 ) 



8 4 40 

+-P(a; 1 )D 1 (t 2 ) + -P 2 (x 1 )D {h) + 2Z>i(t a ) - A^) 



5D (t 2 ) 



+ [ti «-> * 2 ,f!l ^2] 



+ 



Xi + X 2 



C 2 



16D (ti)Do(t2) 



32(D 1 (t 1 )Do(*2) + ^o(ti)Di(t 2 )) 



+ 48A)(*i)A)(*2) 



L M + 32D 1 (t 1 )D 1 (* 2 ) 



+ 24(P 1 (t 1 )Po(t 2 ) + A)(*i)£>i(t 2 )) + (17 + 80C(2))D (ti)D (*2) 



/268 



44 



44 

L M - — (£>i(ti)D (* 2 ) + A)(ti)ZMt 2 )) 



+ 



V 9 



8((2) A)(*i)A)(* 2 ) 



+n f C F 
40 



-D (ti)^o(t 2 ) 



-yA,(*i)A>(k)} 



(B.5) 



If one chooses the renormalization scale hr unequal to the mass factorization scale \x 
one has to add the following term to the expressions in ( |B.1| ) and ( |B.5| ) 



"sOD , tfk A S+V,(1) ( . , n2 2> 

47T /X 2 



where po is the lowest order coefficient in the /3-function given by 

r ll r 2 
Po = y<-A - gW/ 



(B.6) 



(B.7) 



and A^ +y,(1) can be found in (PP[ ) for the MS-scheme and in (|A.13|) for the DIS- 
scheme. 

The coefficient functions for the cross section d 2 a/dQ 2 dy can be very easily derived 
from the above expression by multiplying the coefficient functions in (p.l| ),( |BT5|) and 
( |B.6| ) by the factor x\ + x 2 . 
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Figure captions 

Fig. 1 The ratios R m (y/f) flU) and R {2) (y/f) Q presented in the DIS-scheme. 
Solid line: R {1) (y/f)\ dotted line: R [2) (y/f). Top window : VS = 38.8GeV/c 
(0.125 < y/f < 0.342, E772). Middle window: v 7 ^ = 21.8GeV/c (0.185 < 
y/f < 0.575, E615). Bottom window: VS = 15.4GeV/c (0.26 < y/f < 0.60, 
E537). 

Fig. 2 The same as in fig. [l| but now for the MS-scheme. 

Fig. 3 The ratio dR {l) (y/f,xp) ( p.5|) presented in the DIS-scheme for tc + W — > 
/i+/x~ + "X" at = 15.4GeV/c (E537). Solid line: y/f = 0.25; dotted line: 
y/f = 0.42; dashed line: y/f = 0.60. 

Fig. 4 The ratio dR^ (y/f, xp) ( |3.5|) presented in the DIS-scheme for tc~ + W — > 
fi + fi~ + "X" at \/S = 21.8GeV/c (E615). Solid line: y/f = 0.18; dotted line: 
y/f = 0.42; dashed line: y/f = 0.65. 

Fig. 5 The ratio d R^ (y/f , x p) (|3.5f) presented in the DIS-scheme for p + 2 iJ — > 
+ "X" at v^S = 38.8GeV/c (E772). Solid line: y/f = 0.13; dotted line: 
y/f = 0.23; dashed line: y/f = 0.34. 

Fig. 6 The same as in fig. but now for the MS-scheme. 

Fig. 7 The same as in fig. |4] but now for the MS-scheme. 

Fig. 8 The same as in fig. but now for the MS-scheme. 

Fig. 9 da I dm (|]|) for the reaction p+W -> / u + /i- + "X" at VS = 15.4GeV/c. The 
data are obtained from the E537 experiment |2(|. Dashed line: Born; dotted 
line: 0(a s ) (DIS); solid line: 0(a 2 s ) (DIS); long dashed line: 0(a 2 ) (MS). 

Fig. 10 The same as in fig. |9| but now for the reaction 7r~ + W — > + "X". 

Fig. 11 da/dxp (0) for the reaction p + W -> + "X" at v 7 ^ = 15.4 GeV/c. 

The data are obtained from the E537 experiment [p0[| . Dashed line: Born; 
dotted line: 0(a s ) (DIS); solid line: O(a^) (DIS); long dashed line: 0(a 2 ) 
(MS). 

Fig. 12 The same as in fig. |ll| but now for the reaction ir~ + W — » + "X". 

Fig. 13 da/dy/f = VSda/dm (see (CT) for the reaction 7T + -> fM+p,' + "X" 
at = 21.8GeV/c. The data are obtained from the E615 experiment pT 



Dashed line: Born; dotted line: 0(a s ) (DIS); solid line: 0(a 2 ) (DIS); long 
dashed line: 0(a 2 ) (MS). 

Fig. 14 d 2 a/dy/fdx F QjDj) with 0.185 < y/f < 0.231 for the reaction ir~ + W -> 
+ "X" at VS = 21.8GeV/c. The data are obtained from the E615 
experiment fl2"l| . Dashed line: Born; dotted line: 0(a s ) (DIS); solid line: 0(a 2 ) 
(DIS); long dashed line: 0(a 2 s ) (MS). 
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Fig. 
Fig. 
Fig. 
Fig. 
Fig. 
Fig. 



15 The same as in fig. ITJ| but now for 0.231 < y/r < 0.277. 



16 The same as in fig. [14| but now for 0.277 < \fr < 0.323. 

17 The same as in fig. ITJ| but now for 0.323 < \pr < 0.369. 



18 The same as in fig. [L4| but now for 0.369 < y/r < 0.484. 



19 The same as in fig. |T| but now for 0.484 < y/r < 0.599. 



20 Order a\ corrected K-factor denoted by Ki ( |3.9| ) compared with the ex- 
perimental .ff -factor ( 3.10|) for the reaction ix~ + W — ► + " X" at \/S = 



Fig. 



Fig. 



21.8 GeV/c. The data are obtained from the E615 experiment [2TJ. Dotted line: 
Kx (DIS); solid line: K 2 (DIS); long dashed line: K 2 (MS). 

21 m 3 d 2 a/dmdx F for the reaction p + 2 H -> / u + /i~ + "X" at y/S = 38.8 GeV/c 
and m = 8.15GeV/c 2 . The data are obtained from the E772 experiment p3| . 
Dashed line: Born; dotted line: 0(a s ) (DIS); solid line: 0{a 2 s ) (DIS); long 
dashed line: 0(a 2 s ) (MS). 

22 Parton density dependence of m 3 d 2 a/dmdxp corrected up to order a 2 for 
the reaction p + 2 H -> fi+fi' + at VS = 38.8 GeV/c and m = 8.15 GeV/c 2 . 
The data are obtained from the E772 experiment p3| . Solid line: MRS (SO); 
dotted line: MRS (DO); dashed line: MRS(D-). 
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